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Abstract 

Except for the universe, all quantum systems are open, and ac- 
cording to quantum state diffusion theory, many systems localize to 
wave packets in the neighborhood of phase space points. This is due 
to decoherence from the interaction with the environment, and makes 
the quasiclassical limit of such systems both more realistic and sim- 
pler in many respects than the more familiar quasiclassical limit for 
closed systems. A linearized version of this theory leads to the correct 
classical dynamics in the macroscopic limit, even for nonlinear and 
chaotic systems. We apply the theory to the forced, damped Duffing 
oscillator, comparing the numerical results of the full and linearized 
equations, and argue that this can be used to make explicit calcula- 
tions in the decoherent histories formalism of quantum mechanics. 



QMW Theory Preprint QMW-PH-95-33 

1 



1 Introduction 



One of the chief questions in quantum theory is how classical behavior arises 
in quantum mechanical systems. A naive description of a closed system, 
using the Schrodinger equation, very quickly loses all resemblance to classical 
mechanics. Even a highly localized wave packet can, given sufficient time, 
spread until it occupies all accessible regions of phase space. Such delocahzed 
states are highly non-classical, and quite unlike what we normally observe. 

However, if there is any significant interaction with the environment, such as 
a measurement, the system becomes open, and the Schrodinger equation is 
no longer valid. Macroscopic systems, with the exception of the entire uni- 
verse, are never completely isolated. Even in interstellar space, interactions 
with thermal background radiation and encounters with stray atoms prevent 
macroscopic superpositions from arising [|T|. 

The master equation approach has long been used to treat open quantum 
systems. In it, one averages over the unknown effects of interaction with the 
environment to compute the probabilities of possible output states. This has 
been used in an enormous range of experimental set-ups. 

Modern experiments, however, often enable one to follow the behavior of a 
single quantum system as it evolves, influenced (but not destroyed) by the 
environment, including the measurement apparatus itself. Quantum state 
diffusion (QSD) provides a formalism which can be used to describe such 
situations 0. In it, a system is always in a pure state, evolving according 
to a stochastic nonlinear Langevin-Ito equation. The equation includes both 
collective effects of the environment (such as dissipation and localization) 
and random fluctuations. For most open systems, the classical limit is qual- 
itatively very different than for closed systems. QSD shows this particularly 
clearly, with its dynamical localization producing wave packets which follow 
approximately classical trajectories. In this sense, the classical limit of open 
systems is far simpler than that for closed systems. 

In section 2 we describe the QSD formalism, and compare it with the de- 
coherent histories approach to quantum mechanics ^, ^, ^ |^, ^. It has 
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recently been shown that these two approaches are related, and in fact pro- 
vide complementary views of systems interacting with environments 



In section 3 we examine the quantum dynamics of linear systems and their 
classical limit, and how nonlinear systems can also be treated in this limit. 
We show that this is very different for open systems with localization than for 
closed systems. We briefly outline the importance of localization in efficient 
numerical simulation of QSD, and how it gives QSD remarkable advantages 
over other simulation methods for open systems. This builds on the results 
in 



T0|] , for a particular example. 



In section 4 we examine the problem of quantum dissipative chaos. While 
considerable work has been done on quantum equivalents to classical Hamil- 
tonian chaotic systems, very little has been done for dissipative chaos. In 
part, this is due to the challenging problems involved in studying dissipa- 
tion in quantum systems. Spiller and Ralph |]ll| have shown how quantum 



state diffusion can be used to model dissipative chaotic systems. We exam- 
ine a particular model, the forced, damped Duffing oscillator, to illustrate 
both the classical and quantum characteristics of dissipative chaos, and argue 
the relevance of these results to treatments using decoherent histories. The 
theoretical arguments are backed up by numerical results. 

Finally, in section 5 we re-examine our results, and draw conclusions on the 
usefulness of QSD both as a theoretical and a practical tool in studying the 
quasiclassical limit of open quantum systems. 



2 Quantum dynamics 

2.1 Quantum state diffusion (QSD) 

The master equation is used in the Markov approximation. It can then be 



expressed in Lindhlad form [|12 



z " / 1 1 

P = -J-[H,P] + T.[L,pL] - -L]L,p - -pL]L,), (1) 
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where p is the density operator, H is the Hamiltonian and the m operators 
Lj are the Lindhlad operators, which model the effects of the environment. 
Different choices of these Lj have been used to simulate a wide range of 
physical situations, including measurements. 

The evolution of a single open system is unpredictable. While an ensemble 
of systems may evolve in a perfectly deterministic manner, any single ex- 
perimental run must have a strong stochastic component. There are many 
ways of "unraveling" the master equation into trajectories; the quantum 
state diffusion (QSD) formalism is one of these, with unique symmetry and 
localization properties P, 

In QSD, a state evolves according to the Ito equation 

\d^) = --H\m + T.y^L])L, - -Cl]){L,) - -L]L,p)dt 

m 

+ J2iL, - {L,)md^, , (2) 
i=i 

where the dC,j are m independent complex random differential variables rep- 
resenting a complex Wiener process. They obey relationships 

M{d^,)=M{d^idQ = 0, M{diidi*) = dt5i,, (3) 

where M represents the ensemble average. 

This unraveling of the master equation has many advantages. It is invariant 
under unitary transformations among the Lindblad operators, an invariance 
shared with the master equation, and which is unique to QSD. By taking the 
ensemble mean of the projector \ip){ip\, one reproduces the master equation 
above: M {\%p) = p. This gives a very simple interpretation for the rela- 
tionship between the trajectories li'it)) and the density operator p: p(t) is 
the ensemble of all possible states over the course of many runs of an 

experiment, the \ip)^s will tend to reproduce the density operator p, while in 
a single run, the system is in a state {ip) (which depends on the choice of un- 
raveling). Thus, expectation values {0)p calculated with a density operator 
p would include both the quantum expectation value (O)^ and the ensemble 
mean over all the \ip)s. 
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Numerically there is often a large advantage in representing a system by a 
state rather than a density operator. If it requires basis states {2N real 
numbers) to represent the wavefunction, the density operator will require N"^ 
real numbers. Thus, even though many runs may be required to produce an 
accurate answer, the factor of order N saved on each individual run can still 
result in large savings of computer time and memory. This is particularly 
true if one takes advantage of QSD's localization properties fTO . 



2.2 Decoherent histories 



An extremely promising interpretation of quantum mechanics is the decoher- 
ent (or consistent) histories formalism, developed by Griffiths, Omnes, and 
Gell-Mann and Hartle [^, ^, |], |^, ^ . In this formalism, quantum systems are 
described by sets of possible histories {a}. These histories ascribe certain 
definite values, or ranges of values, to selected variables at successive points 
in time. Along with the histories, there is a decoherence functional D[a, a'] 
on pairs of histories. For a set of histories {a} to consistently describe the 
quantum system, they must satisfy the decoherence criterion 

D[a,a'] = Saa'Pa, (4) 

where pa is the probability for the history a actually to occur. If the decoher- 
ence criterion is met, then these histories obey the usual classical probability 
sum rules; one can assert that one or another history actually occurred, with- 
out the kind of interference effects seen in the two-slit experiment. In the 
case of Schrodinger's cat, one history would include a live cat and another a 
dead cat, without reference to observers or measurements. 

In standard non-relativistic quantum mechanics, the usual choice of histories 
consists of an exhaustive set of orthogonal projection operators {P^. (ti)} at 
each time of interest tj. 

T.KM = l ^:.('^.)K^('^.) = ^.,.^n(tO- (5) 

These projections represent different alternatives at the time tj. A single 
history a represents a choice of one Vl^.{ti) at each time tj, represented by 
the string of projection choices a^. 
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The decoherence functional can be written 



D[a,a'] = Tv\viStn) ■ ■ ■'Pl{h)pr\{t,) ■ ■ -P^, (tn)}, (6) 
and a set of histories must satisfy (^. 

One might be tempted to choose as a set of histories the precise values of 
some complete set of variables at each point in time. For example, for a single 
particle one might have its histories be its different possible trajectories x{t). 
Histories of this sort almost never decohere. In general, a fair degree of 
coarse-graining is required to even approximately satisfy (^. Since finding 
exactly decoherent sets of histories is difficult, the usual course is to find a 
set of histories which satisfies the decoherence criterion to a finite level of 
accuracy fl^. It is conjectured that there is an exactly decoherent set of 



histories in a neighborhood of every such approximately decoherent set [15 



One simple (and therefore popular) type of coarse-graining is to separate the 



system into a distinguished subsystem and an environment ^, 16|. Tracing 
over the environment then produces a set of coarse-grained histories in terms 
of the "distinguished variables" alone. It has been shown that under a wide 
variety of circumstances, histories coarse-grained in this way approximately 
decohere, with increasing precision as the distinguished system approaches 
the classical limit [Q. 

This type of coarse-graining brings us close to the situation of interest for 
quantum state diffusion. It has recently been shown that for a large class 
of open quantum systems there is a close connection between quantum state 
diffusion and decoherent histories 0. The different possible QSD trajectories 
correspond to alternative histories; these histories approximately decohere, 
and occur with the same probabilities as in the decoherence formalism. In 
this context one can think of QSD as a dynamical form of the decoherent 
histories formalism, or decoherent histories as a global view of QSD. 



3 Open quantum systems and the quasiclas- 
sical limit 
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3.1 Linear quantum systems 

An important class of systems which can be treated with the QSD formahsm 
are linear systems, with a Hamiltonian function quadratic in Q and P, and 
Lindblad operators given by a (not necessarily Hermitian) combination CqQ + 
CpP. These systems have been treated at length by a number of researchers 

MM- 

Halliwell and Zoupas [|19] have shown that these linear systems tend to lo- 
calize with time onto Gaussian wave packets in phase space. These wave 
packets follow paths in phase space which are close to paths defined by the 
classical (or quasiclassical) equations of motion. 

From the QSD equation (H) the Ito equation for the expectation value of an 
observable is 

d{d) = (#|6|^) + {ip\d\diij) + {diij\d\dij) 

= -i{[d, H])dt + {DOL - i{L^ L,d})dt + {{OL) - {d){L))d^ + 

((6Lt)-(6)(Lt))de- (7) 
If the Hamiltonian is of the usual form 

H = Py2m + V{Q), 

and the Lindblad operator L is linear then we can readily evaluate this ex- 
pression for simple operators O. 

Halliwell and Zoupas have shown that for linear systems one can, to a very 
good approximation, reduce the QSD equation to a set of coupled equations 
for (Q), (P), (Aq\ (Ap\ and {AQAP + APAQ). So long as the wave 
packet remains localized, so that the higher moments can be neglected, the 
QSD evolution can be reduced to these five equations. The equations for {Q) 
and (P) correspond to the equations of motion for a classical path in phase 
space, modified by the presence of noise and small diffusive corrections which 
depend on the higher moments. 

This is similar to the Ehrenfest theorem for the Schrodinger equation, in 
which it can be shown that wave packets localized in phase space follow 
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approximately classical paths. Unlike the Ehrenfest theorem, however, in 
QSD this localization usually becomes better and better with time, rather 
than the reverse |]TB|, |TU|. 



3.2 Nonlinear quantum systems 



While the equations of Halliwell and Zoupas were derived for linear systems, 
they provide a good approximation for the classical limit even of nonlinear 
systems. In this classical limit, the width of a localized wave packet is small 
compared to the spatial variance of the potential (and its higher derivatives). 
Thus, to the wave packet the potential seems locally quadratic, though slowly 
time- varying, and we would expect many of these localization results to apply. 

In the classical limit, therefore, one can treat nonlinear problems by lineariz- 
ing about the quasiclassical path, and solving for the moments of the wave 
packet. This is quite similar to the classical technique of linearizing about a 
known orbit to examine the behavior of nearby solutions. 

Halliwell and Zoupas [0 have shown that solutions tend in the mean towards 
localized Gaussian wave packets when the third-derivative and higher terms 
in the potential are negligible, i.e., when 

{v'm ^ v'm\ (8) 

{QV'm - {Q){V'{Q)) ^ (AqV((Q)), (9) 

and 

{PV\Q) + V\Q)P) - 2{P) {V'iQ)) ^ {AQAP + AP AQ)V" {{Q)) . (10) 



Numerical investigation |g, |T0[ as well as theoretical arguments [|13l indicate 
that localization occurs in almost all systems, whether linear or not. If one 
thinks of QSD in terms of decoherence this is not surprising. Many people 
have shown that decoherence through interaction with an environment is 
extremely rapid. This corresponds to superpositions quickly reducing to 
highly localized states. In most cases, this localizing effect of the environment 
predominates over dispersive effects due to quantum spreading. 
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3.3 Localization and the moving basis 



One practical consequence of localization is that it is possible, by choosing an 
appropriate basis, to reduce significantly the number of basis states needed 



to represent the wave packet |10|. If a wave packet is localized about a point 
{q,p) in phase space, it might require a great many of the usual number 
states \n) to represent it. If instead we represent the wave packet in terms 
of excited coherent states \q,p,n) = D{q,p)\n), relatively few basis states are 
required, with obvious savings in terms of computer storage and computation 
time. Here, D[q,p) is the coherent state displacement operator, 

D{q,p)=exp^i^pQ-qpy (11) 

which displaces the excited state \n) to the excited coherent state \q,p, n) cen- 
tered on This representation of the state by classical and quan- 
tum ra) is called the moving basis or mzxerf representation, or MQSD, 
and has been discussed extensively elsewhere [p^, . 



4 The Duffing oscillator 



4.1 Quantizing the forced damped Duffing oscillator 

A good example of a nonlinear system is the forced, damped Duffing oscilla- 
tor. This has a classical equation of motion 

d X „ _ doc Q . > / ^ _ \ 

— + 2r— + x''-x=g cos(t), (12) 



and for some choices of F and g is chaotic pT 



Because the equation of motion includes explicit time-dependence, the so- 
lutions lie in a three-dimensional phase space x,p,t. It is helpful to con- 
sider a discrete surface of section of this system. Let (xcPo) be the initial 
point of the forced, damped Duffing oscillator at time to = 0- Then we 
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can define a constant phase map in the x-p plane by the sequence of points 
{xn,Pn) = {x(tn),p(tn)) at timcs tn = 27m. Figure 1 illustrates this in the 
chaotic regime, where we can clearly see from the surface of section the fractal 
structure of the strange attractor. 

Quantizing the Duffing oscillator is straightforward using the QSD formalism. 
The Hamiltonian operator is 

HiQ, P, t) = P^/2m + QV4 - + 9 cos{t)Q + Vf{Qp + PQ), (13) 

and the damping is represented by a Lindblad operator 

L = 2Vfa = V2r{Q + tP), (14) 

where we have assumed h = 1. The last term in the Hamiltonian is an 
ansatz, added to give the correct equations of motion in the classical limit; 
it is necessary due to the simplistic model of the dissipative environment. 

This system is far from classical. To go to the classical limit, we introduce a 
scaling factor (3, 

HfsiQ, P, t) = Py2m+p^Qy4:-Qy2+{g/p) cosit)Q+Vf{Qp+pQ). (15) 

As we reduce (3, the scale of the problem (compared to h) increases by 1//3 in 
X and p, without altering the classical dynamics. Thus, P ^ is the classical 
limit of this system. Classical behavior should emerge from the system in 
this limit. This is supported by the numerical calculations. (See section 4.3.) 



4.2 Application to decoherent histories 

The forced, damped Duffing oscillator has been treated elsewhere using the 



decoherent histories formalism [^, |2^, While this approach has great 
theoretical power, it does not lend itself readily to numerical solution, which 
is generally necessary for this kind of nonlinear problem. Fortunately, as it 
has been shown that QSD is very closely related to a certain kind of coarse- 
grained decoherent history 0, one can use the numerically efficient QSD 
method to perform a quantum Monte Carlo calculation for open systems of 
this sort. 
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It has been argued |T^, ^ |T^, ^ that solutions of these open systems should 
peak about the quasiclassical equations of motion, i.e., they should behave in 
the classical limit like a classical problem with noise. This is amply demon- 
strated by the numerical results in section 4.3. 

By choosing an appropriate set of decoherent histories one can find a quantum 
equivalent to the Poincare surface of section defined by (|12D above. Consider 
the times ti = 27ci, with a set of projection operators Vq-^pXti) onto local cells 
of phase space defined at each time, these cells centered at {qi,Pi)- Although 
there are no true orthogonal projections onto cells of phase space, there are 
a number of approximate projections; the simplest are the coherent state 
projectors Vqp = \q,p){q,p\, and there are many others P, |2^, ^ ^ . 
We can write down the decoherence functional 

D[{q^,p^}, = Tr|p,„,p„(t„) • ■ ■K,p,ih)pVq'^,p>ih) ■ ■ ■ P,; 

(16) 

where p = is the initial density operator, for simplicity taken to be 

pure. 

It has been shown that in systems which interact strongly with an environ- 
ment, the off-diagonal terms of the decoherence functional decay very rapidly. 
So long as the time between projections tj+i — tj is large compared to the 
decoherence time to, and the phase space cells are large compared to h, the 
histories described by ( |16|) will decohere to a good level of precision. In the 
quasiclassical limit, the probabilities D[{qi,pi}, {qi,pi}] should be sharply 
peaked about the classical solution. 

As shown by Diosi, Gisin, Halliwell and Percival |0, quantum state diffusion 
trajectories will give the same probabilities as the decoherence functional 
for histories which project onto the same quantities (in this case coherent 
states in phase space) and are at least as coarse-grained as the level of QSD 
localization. By following a QSD trajectory, therefore, and looking at the 
expectation values of {Q) and (P) at the times ti, one has picked out a single 
decoherent history with the correct probability. 
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4.3 Numerical example 



We solved the forced, damped Duffing oscillator numerically by integrating 
the QSD equation, both in full (using MQSD) and in the linearized approx- 
imation of Halliwell and Zoupas, using a C++ quantum simulation package 
we have developed [|IU], All of these results were generated on standard 



486 PC's running LINUX. 

The classical unsealed problem is bounded within a small region of phase 
space. (See figure 1.) Clearly, the quantized version of this problem (with 
h = 1) should be far from the classical limit. One would expect to observe 
little trace of the classical fractal structure. The expectation values of Q and 
P should be dominated by noise. 

Examining figure 2a, we see that the results of the numerical calculation 
match our expectations very well. The expectation values appear randomly 
distributed; they are dominated almost completely by the stochastic terms 
of the equation. 

Approaching the classical limit /5 — 0, more and more of the classical struc- 
ture of the attractor appears. At ffist the broad outlines of the attractor 
are formed, then increasing levels of substructure. (Figures 2b-2d.) The full 
fractal structure of a strange attractor is only attainable in an unphysical 
purely classical limit. In an actual physical system, the uncertainty principle 
provides a lower cutoff to the scale-invariance of the strange attractor. It is 
unphysical to discuss regions of phase space with areas smaller than Planck's 
constant. 

In fact, the actual lower cutoff is even larger than that, due to the interaction 
with the environment. Classically, one can eliminate noise from the environ- 
ment by going to the limit of zero temperature. In quantum systems, though, 
there is noise even in this limit - the equivalent of zero point noise for the 
degrees of freedom of the reservoir. The existence of this noise is closely re- 
lated to the fact that in general, coarse-grainings much coarser than would be 
naively suggested by the uncertainty principle are required for decoherence 
P^ ; in general, phase space cells large compared to h are required. 
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4.4 The linearized theory 



The linearized equation are expected to remain valid as long as the potential 
is approximately locally quadratic, i.e., as /3 — > 0. This approximation takes 
advantage of QSD's unique localization properties. We can solve these equa- 
tions, suitably modified by the presence of the ansatz in (|T^), for different 
values of /?, and compare them to results generated with the full QSD code. 
The results are displayed in figures 3a-3d. 

We see that at /5 = 1 the linearized equations are completely wrong (figure 
3a). As we reduce (3, however, the two methods come into better and better 
agreement, until at = 0.01 the agreement is nearly perfect (figures 3b-3d). 

The linearized equations are thus only valid when one goes very far toward 
the quasiclassical limit. In this regime, however, the computational advan- 
tage over the full equations is enormous; for /3 = 0.01 it is 10^. In the 
intermediate regime of (3 ~ 0.1, the linearized approximation is no longer 
valid, but localization can still provide a tremendous advantage. Using the 
MQSD algorithm, one can represent the system by 10-20 basis states, as 
opposed to the ~ 1000 basis states which would be needed by QSD without 
the moving basis algorithm and more than 10^ real numbers which would be 
required to solve the master equation. 



5 Conclusions 



Many criteria have been suggested to define chaos in quantum systems. These 
differ in their emphasis and domain of application. But nearly all are based 
on the premise that, as one approaches the classical realm, quantum systems 
should behave more and more like classical chaotic systems. This limit is 
difficult in the closed systems which are usually studied. By contrast, for 
open systems the classical limit arises naturally. 

Quantum state diffusion provides a vivid picture of how this crossover from 
quantum to classical occurs. As a quantum system becomes macroscopic, it 
increasingly resembles a localized wave packet, following an approximately 
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classical path in phase space. At short length scales, its evolution is domi- 
nated by quantum mechanical uncertainties, due to quantum spreading and 
the random influence of the environment. As one approaches the classical 
limit, these stochastic influences become less important, being replaced by 
the uncertainties of classical chaos. 

This is clearly illustrated by the numerical results of QSD simulations. These 
demonstrate the value of QSD as a practical algorithm for solving the mas- 
ter equation. By making use of QSD's unique localization properties, the 
equations can be solved numerically with far greater efficiency than would 
be possible in solving the master equation directly [ p!0[] . 

Going to the classical limit this advantage is even more pronounced; one may 
replace the full QSD equations by a linearized approximation, requiring the 



solution of only 5 coupled ordinary differential equations |T9[. This approx- 
imation provides not only a powerful numerical tool, but also demonstrates 
how quasiclassical equations of motion rise from the underlying quantum 
theory. 

QSD also provides a practical method of calculating results in the decoher- 
ent histories formalism of quantum mechanics. While decoherent histories 
provides a powerful interpretive tool, it does not necessarily lend itself to 
practical calculations; it is often difficult to solve for the probabilities of 
individual histories. For open systems, QSD generates such histories au- 
tomatically, with their correct probabilities. The arguments of Gell-Mann 
and Hartle, that probabilities will peak about quasiclassical histories in the 
classical limit, are supported by the QSD analysis of this limit, and by the 
numerical results presented in this paper. 
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Figure 1. The constant phase surface of section for the classical forced, 
damped Duffing oscillator in the chaotic regime, F = 0.125, g = 0.3. 



Figure 2. The constant phase surface of section for a single QSD trajectory 
of the quantum forced, damped Duffing oscillator in the chaotic regime, F = 
0.125, g = 0.3, for four scalings: a) /3 = 1.0, h) /3 = 0.25, c) /3 = 0.1, d) 
/3 = 0.01. 



Figure 3. The constant phase surface of section for a single linearized QSD 
trajectory of the quantum forced, damped Duffing oscillator in the chaotic 
regime, F = 0.125, g = 0.3, for four scalings: a) /3 = 1.0, b) /5 = 0.25, c) 
/3 = 0.1, d) /3 = 0.01. 
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Figure 2a. QSD solution. 
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Figure 2c. QSD solution. 
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Figure 2d. QSD solution. 
1 



60 - 



40 - 



20 - 



- 



-20 - 



-40 - 



-60 



-150 



o 



o 



o 
o 



o 

$ 

o 



o o 



o mo 



o 



o 



o 
o o 



o 



o 



o 



o 



o 



o 



o 



oo 



o 



o \ 
o o 



o 



-100 



-50 




X 



50 



o 



o 

o 



o 



o 

o 
o 
o 



o 



o 



100 



150 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 

X 



Figure 3d. Linearized QSD solution. 
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